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ABSTRACT 

Bath ythermcgraph data acquired from rhe research vessel 
USNS SILAS BENT along a meridional tracic in the Northeast 
Pacific during September 1977 were statistically analyzed to 
determine possible a ssociaticns between the subsurface ther- 
mal structure and sea surface temperature. Strcngly corre- 
lated variables (thermocline gradients, mixed layer depth, 
and locations cf the seascnal and main thermoclines) within 
the vertical temperature profile were used in linear regres- 
sion methods tc form empirical relationships. The generated 
equations then are utilized to define the subsurface thermal 
structure from only an input of sea surface temperature. 
Comparison tests with temporally and spatially removed BT 
data ware conducted with results indicating successful 
application within a water mass domain with uniformly chang- 
ing character ist ics. 
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INTROD UCTION 



A. IMPORTANCE OF DETERMINING OCEAN THERMAL STRUCTURE 

An important physical property and product of various 
dynamical processes occurring in the oceans is the thermal 
structure. Although it is derived from only a few mecha- 
nisms such as air-sea thermal energy exchange, forced and 
free convective mixing and horizontal ocean thermal advec- 
tion, the interaction of these processes along with other 
hydrodynamic motions can become immensely complicated. 
Since the temperature profile has important implications for 
climatic weather prediction, commercial fisheries and mili- 
tary acoustic surveillance, to mention a few, the ability 
to predict this structure is of major importance. 

Climatological atlases present an averaged profile for a 
certain period and location. However, this results in a 
generalized thermal structure which is not accurate enough 
to accomplish the precise needs of underwater sound propaga- 
tion forecasting. Recent developments in numerical modeling 
of the upper-ocean temperatures have produced vast improve- 
ments in the analysis and prediction of mixed layers and 
other synoptic oceanic features (Clancy et al. , 1S81). It 

is inevitable that this method will reach the sophistication 
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cf today's mstscrclogica 1 forecasting, but wirh the complex- 
ity of the physical processes involved and with present data 
sparsity ir. many ocean areas, many inaccuracies still exist. 

Ihe tathythermo graph (BT) provides an accurate 
temperature profile at a point location in the upper ocean. 
Since the vertical thermal structure over a large area pro- 
vides the main factor which affects the propagation charac- 
teristics cf underwater acoustics, the need for many 
temperature profiles in Navy asw operations is obvious. The 
expense involved in ship and aircraft deployment of BT's for 
an ocean area of interest often prohibits the use of this 
most accurate means of obtaining a detailed thermal 
description. 

B. SATELLITE DE TEH M INA TI C N OF SUBSURFACE THERMAL STRUCTURE 

The study cf the oceans by satellites nas become a major 
arena for scientific scrutiny and investigation. Already it 
has given added insight into the numerous processes that 
govern oceanic behavior. The Gulf Stream eddies and fronts, 
the coastal upwelling off the Somali and Arabian coasts, and 
other applications in biological, geophysical and dynamical 
fields are becoming well known. However, in spits of great 
progress, satellite remote sensing techniques for observing 
oceanic processes are far from being fully realized. 
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The fundamental problem which continuously arises, 
besides atmospheric effects which will not be discussed, is 
that satellite ctservaticns are made at the surface alcne 
and that all ccnclusions about subsurface processes are 
obtained from inference rather than by direct measurement. 
Consequently, the applications cf satellites for oceanogra- 
phy are limited to those phenomena that generate surface 
signatures detectable by electromagnetic waves. Once these 
surface manifestations are observed, it is necessary to have 
knowledge cf the physical processes involved to interpret 
their meaning. Currently there are few methods available 
where a qualitative "feel" for subsurface features can be 
directly expressed as empirical relationships from remotely 
sensed measurements at the sea surface. If a linkage is 
found to exist between the air-sea interface and the thermal 
structure belcw, a satellite would then be the perfect plat- 
form to collect sea surface data rapidly over all the 
world's oceans. This surface information could then be used 
to infer the subsurface vertical temperature structure. 

C. EUHPOSE AKE GCALS OF THESIS 

The primary purpose fcr the cngoing research into the 
prediction of subsurface thermal structure is tc understand 
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the processes and relationships than govern the ocean's 
behavior. Holes that are performed by ships and aircraft 
(via expendable bathythermographs) in determining the 
temperature profiles could be played better by satellites, 
provided a "link" is found to read rhe below-s urf ace events 
from the sea surface signatures. Therefore, if strong cor- 
relations and relationships could be ascertained between 
satellite- derived sea surface temperatures and vertical 
temperature profiles, then the areas of acousrical oceanog- 
raphy and naval tactioal applications would benefit greatly. 

The goal of this thesis is to determine the possible 
correlations that exist between the subsurface thermal 
structure and the temperature at the surface. Then, from 
these relationships and use of regression formulas, 
predicted vertical temperature profiles could be computed 
from an input of sea surface temperature alone. 
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II 



. CCEANOGEAPHIC CHARACTERISTICS OF THE 
121221131 FTOIH'iiiii 

A. KIND SISTEKS AND CURRENTS 

The Northeast Pacific has been studied extensively for 
many years and it has revealed a complex arrangement of 
environmental conditions. Severe storms, high winds, varied 
water masses, numerous currents and unique temperature and 
salinity structures are a normal occurrence for the region. 

In the winter the North Pacific is usually under the 
influence of the Aleutian Low in the atmosphere while for 
the summer months a high pressure system dominates (Uda, 
1963) . The distribution of the mean surface winds are con- 
trolled by these two systems. At the location of the now 
retired Ocean Weather Station ”P” (OWS-P, positioned at 50N, 
145W in the Gulf of Alaska) the prevailing winds are from 
the southwest (Tafcata, 1S65) . From the influence of these 
wind patterns, the surface waters flow eastward across the 
Pacific Ocean in the vicinity of Station "P”. This flow 
consists of the West Wind Drift and Subarctic Current which 
then diverge to form the southward flowing California Cur- 
rent and the northward flowing Alaskan Current (Tabata, 
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1978) 



Fig 



1 shows the location of Ocean Weather Station 



II p II 



and nhe 



major surface 



currents of the 



North Pacific. 




Figure 



Location of OWS-P and North Pacific Ocean 
Currents (From lully, I960). 



Surface 
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B. WATER HASSES AND FRONIS 

The Ncrtheast Pacific can be divided into three major 
regions based upon their water characteris- ics . Near the 
viciniry cf CWS-P and northward is the Subarctic Pacific 
Water Mass in which the upper 100-150 m layer is composed of 
cold and low-salinity water. To the south warm and highly 
saline water exists in the upper layers of the Subtropic 
Water Mass. In between and along the coast lies a region in 
which characteristics frcm both water masses are combined. 
This is known as the Transition Zone where there is a marked 
gradation of oceanographic properties from one region to 
another, with the boundaries separating these water masses 
having complicated therrachaline structures (Roden, 1970). 
The location cf the Subarctic Front between the two water 
masses remains fairly constant, from 40N to 45N, and is 
largely dependent on the prevailing winds (Roden, 1975). 

1. SUBARCTIC WATER MASS CHARAC TER ISTICS 

The main feature of the vertical structure cf 
oceanic prcperties is the existence of an upper zone, a 
halccline and a lower zone (Uda, 1963). The temperature, 
salinity and density profiles of these three zones for sum- 
mer and wintsr are depicted in Fig. 2. 
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Figure 2. Characterisr ic Temperature and Salinity Profiles 
(From Tull y, 1 S64) . 
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The upper zone. 



which extends from the surface tc a 



depth of about ICO m, generally has a low salinity level of 
less than 33 g/kg with the lowest values being closest tc 
the coastal recicns in the late spring <Uda, 1963) . In the 
winter, the water in this zone is thoroughly mixed with all 
physical properties becoming homcgeneous. During the sum- 
mer, due tc lighter surface winds, a shallow ischaline and 
isothermal layer exists from 10 to 30 m below the sea sur- 
face. Below this mixed layer is a strong negative thermo- 
cline with teiperature decreasing as much as 8C in 20 at; in 
the same zone a small halccline extends above another near- 
ischaline layer to the bottom of the upper zone (Tabata, 
1965) . 

Belcw the upper zone lies the main halocline which 
is permanent, as the salinity increases by as much as 1 g/kg 
from 100 tc 200 m (Tabata, 1960) . The temperature generally 
decreases with depth in this region, but at times a 
temperature inversion occurs tc a depth of about 150 m. 
This positive temperature gradient is in a stable layer 
(density increases with depth) since it is located within 
the main halocline (Uda, 1963). 
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In the lower zone, telow the halocline, the proper- 
ties of the water change gradually with depth. Saliniry 
increases slightly by about 0.6 g/kg in 1000 m while the 
temperature decreases to a value of approximately 2.8C at a 
depth of 1000 m (Tabata, 1865). 

The variations of salinity are most pronounced in 
the upper zone and halocline and are less in the lower zone. 
The annual salinity variations in the upper zone are charac- 
terized by a decrease in the upper 30 m in the late summer 
or early fall with a naxiirum in late winter or early spring. 
These seasonal fluctuations result primarily from the evapc- 
raticn and precipitation imbalances common to the region. 
Annual variaticn is similar at great depths, but lags the 
surface layers in time (Tatata, 1965) . 

At the surface, several features are observed in the 
salinity pattern. Basically, the trend of the surface 
isohalines appears to be east-west throughout the year with 
major changes occurring along the coast. Also, a permanent 
salinity maximum cf about 23 g/kg exists near the center of 
the Gulf cf Alaska, south of Kodiak Island. The salinity 
decreases southward to a minimum of about 32.6 g/kg in the 
vicinity of SON and then increases further south into the 
Subtropical Water Mass (Oda, 1963). 
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Th€ ’’three-zcne" structure of the water colum:: is 
characteristic of the deep ocean region of the Subarctic 
Domain. However, it does not always exist in near-coastal 
waters as current systems, coastal winds and fresh water 
runoff from land can modify the profiles. As shown in Figs. 
3 and 4, this coastal regime primarily alters the upper 
zone. Along the western United States this domain exists 
out tc about 130W. The lower zone also becomes affected by 

the proximity to land of the California Undercurrent 
(Dcdimead et al. , 1962). 

2. SUETBCEIC WAT EH MA^ CHARACTERISTICS 

To the south of the Subarctic region is the Sub- 
tropic Water Mass which is basically warmer and more saline 
than the Subarctic Water. In the upper layer the high val- 
ues cf salinity (greater than 34 g/kg) due to higher evapo- 
ration decrease to a minimum at 200 to 600 m and are typical 
of the values throughout the year. Sea surface temperatures 
remain fairly constant with an isothermal upper zone averag- 
ing about 20 m in thickness below the sea surface, due tc 
the prevailing winds. Eelow this zone exists a negative 
permanent thermocline reaching dcwn to a depth of 500-600 ra. 
Use, there is no characteristic halocline or major 



25 




Figure 3. Upper Zone Domains (From Dodimeai ei 



1962) . 
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Figure 4. Lower Zone Domains (From Dodimead ^ al . « 



1 962) . 
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temperature irxversicn as there was in the warer mass to the 
north (Tatata, 1S65). The deep wa-er structure is basically 
continuous and similar from the Subarctic to Subtropic 
regicns (Dda, 1963) . 

3. TRANSIIICN ZONE 

The Transition Zone, located in the central and 
eastern North Eacific between 45N and 35N , separates Subarc- 
tic Water from Subtropic Water. It consists of the mixed 
water originating in the warm, saline Karcshio and the cold, 
low salinity Cyashio currents as they move eastward in the 
West Wind Drift (Uda, 1963). 

The width of this zonal belt, which is about 2 tc 4 
degrees of latitude, and the strength of its boundaries, are 
primarily derermined by the wind stress at nhe sea surface. 
The ncrrhern boundary is characterized by many temperature 
inversions and by the gradual disappearance of -he Subarctic 
halocline. At the southern boundary the vertical structure 
becomes basically isohaline with a strong negative thermo- 
cline in the upper layer (Scden, 1970). 

The Transition Zone is primarily a property of the 
upper ocean, as its char acter ist ics fade below a few hundred 
meters. At the surface this region has a distinctive 
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meridional (N-S) temperature gradient of approximately 1.0C 
per 100km as it stretches eastward across the Pacific Ocean. 
The horizcncal surface salinity pattern exhibits similar 
structure alcng a latitude band throughout the year with 
values increasing as much as 0.8 g/kg meridionally (Uda, 
1963) . 

The upper layers are often characterized by 
temperature inversions which can occur because of the exis- 
tence of the associated halocline. These inversions 
decrease in magnitude and slowly increase in depth toward 
lower latitudes. 

The vertical distribution of properties at interme- 
diate depths shows that the salinity becomes effectively 
isohaline, with a slight minimum developing below 300 m at 
the southern boundary (Uda, 1963). 

STRUCTURE AND THE SUBA^TIC FRONT 

The density of sea water is calculated primarily 
from its temperature and salinity at depth. In the Subarc- 
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temperature and salinity are well mixed. In summer, a minor 
pycnccline corresponding to the shallow rhermocline is 
observed. However, a large and permanent pycnocline occurs 
along with the main halocline. In the lower zone, the den- 
sity increases gradually with depth (Uda, 1963). 

In Subtropic Water, the density is controlled by the 
temperature. Hare the warm surface layer of about 10 to 
30 m is isopycnal from the mixing of the prevailing winds. 
A permanent pycnccline exists at lower depths due to the 
thermccline in the region (Oda, 1963) . 

The Transition Region has a complex density struc- 
ture as numerous minor density fronts occur at depth due to 
the temperature and salinity inversions. 

The Subarctic Front in the eastern North Pacific is 
a narrow, meandering band located in the latitude range of 
40N to 45N. This front has a unique structure in its lack 
of a density front in the upper 100 m. The strong horizon- 
tal temperature and salinity gradients are in almost com- 
plete balance. Below this upper layer a density front of 
moderate intensity slopes slightly to the south (Roden, 
1975) . 



30 



The mixed layer deprh can szrongly change across rhe 
Subarctic Frcrt. Tc the norrh, rhe mixed layer extends to 
the top of the halocline at about 100 m. To the south, 
there is nc halccline and, during the stormy, winter season, 
the mixed layer depth can extend to almost 300 m. In the 
summer months the seasonal thermocline controls the depth of 
the mixed layer and there is little difference across the 
Subarctic Front (Eoden, 1S75). 
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III. THERMAL STRUCTURE 



A. THERMAL STRUCTURE IN THE NORTHEAST PACIFIC OCEAN 

The ifiteracticn between the sea and atmosphere is the 
dominant process determining temperature structure in the 
upper layers. Below, temperature is a fairly consrant prop- 
erty with the structure being determined by internal pro- 
cesses, such as advection and mixing. 

Two heating cycles occur in the ocean; diurnal and sea- 
sonal. In the daily cycle rhe surface layer is heated dur- 
ing the daylight hours and cooled during the nighr. In the 
higher latitudes during the heating season, i.e., the middle 
six months of the year, a daily net heat gain occurs (Uda, 
1963). Fig. 5 illustrates the annual cycle of heating and 
cooling and its relation to the thermal structure. In 
March, at the end of the cooling season, the waters are 
nearly isctheiial to the top of the halocline located at 
about 100 m in depth. With the beginning of the heating 
season in April, the heat is accumulated as small 
temperature inversions are gradually mixed downward by light 
winds. During subsequent short periods of high winds these 
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transi'riits are guiclcly eroded and a shallow isothermal sur- 
face layer results. 

Eelow the surface layer, a sharp negative ther mocline is 
formed between 30 and 60 m in a temperarure transition layer 
adjacent to the deeper cold waters. Prom the end of the 
heating season until the waters are again nearly isothermal 
in March, the thermocline deepens to about 100 m. 

In this ocean region, below the rhe rmocline there is 
usually a cooler layer nearly 120 m in depoh which is sra- 
ble, since it occurs within the main halocline. Below the 
cooler layer there may be an increase in temperature with 
depth, an inversion. Temperature inversions can be classi- 
fied into three types according to their method of formation 
and behavior. The first type, which is the most common in 
the eastern North Pacific, is a result of the seasonal cycle 
of heating and cooling. The upper layers in an extremely 
cold winter can become colder than the waters in the halo- 
cline below and a positive temperature gradient is formed. 
With the advent of the following heating season the 
temperature of the surface waters increase causing a 
(negative) thermocline tc form above the positive thermal 
gradient. The temperature minimum thus generally is located 
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Figure 5. Growth and Decay of the Thermociine at Ocean 
Szaticn ”P" (Ficm Tally and Giovando, 1963). 
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just above the halocline at 120 ra and changes very little 
during the course of the summer heating period. During the 
following winter season the upper thermocline decays and the 
temperature inversion may be erased (Uda, 1963). 

another process which can causa temperature inversions 
is the formation of sea-ice in the western Pacific. Tha 
remaining cold (-1C) , high salinity surface layer sinks to a 
level of stability in the existing halocline and is carried 
eastward in the West Wind Drift where it completely mixes 
with other waters before reaching 140W, Depending cn the 
amount of heating of the surface layers during the transit, 
the temperature inversion may or may not persist throughout 
the year (Eoden, 1964). 

a third process results in a subsurface temperature max- 
imum when the warm, more saline Subtropical Water intrudes 
under the colder Subarctic Waters. The advecting water then 
is slightly cooled and sinks; it remains warmer than the 
surroundings, however, and forms the temperature maximum at 
depth. Its upper boundary is generally limited by the sta- 
bility criterion of the main halocline (Uda, 1963) . 

As the cooling period begins in lata September, the sur- 
face waters are cooled and mixed downward by convection as 
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wall as cy the strong winds of the fall and winter months. 
The seasonal thsrmocline sinks and narrows unt.il it reaches 
the halocline. By late winter the upper zona has reached 
its minimuB reicperature of the year. 

At depths telow 500 m temperature variations are small, 
seldom exceeding standard deviations of 0.04C (Tabata, 
1960). Fig. 6 shows the typical temperature structure at 
Ocean station "F" for both summer and winter months. 

B. INTERNAL WAVE EFFECTS CN THERMAL STRUCTURE 

In many vertical temperature structure analyses, the 
depths of the isotherms fluctuate considerably. Periodicity 
and amplitude of the fluctuations appear non-systemat ic , but 
result from the superposition of a large number of internal 
wave trains. The limiting frequencies range from the 
Br unt-Vaisala to inertial; in the Subarctic Region, periods 
from 5 minutes to at least half a day have been observed. 
The vertical displacement of isotherms is of the order of 
5 m at the top of the thermccline and perhaps as large as 
25 m at the level of the principal halocline (Tully and 
Giovando, 1963) . 
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Figur<= 6. 



Tycical Tempera-urs Structures at Ocean Station 
"pn (From Tully and Giovando, 1963). 
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C. HC2IZ0NTAI SURFACE THEEHAL STRUCTURE 

Near the center of the Gulf of Alaska is located an area 
of temperature ffiniouii that coincides with the dynamic cen- 
ter of the Alaskan Gyre. In summer it roughly overlies the 
center of the gyre as determined by the salinity maximum, 
but in winter the two are separated (Uda, 1963). 

South from the cold core dome to approximately 35N the 
thermal field is oriented zonally except near the coastal 
regions. This configuration continues year round with a 
seasonal range of 7.0C throughout the Transition Region 
(Uda, 1963). The horizontal thermal gradient in the north- 
south direction is strongest in the Transition- Zone where 
gradients of IC/IOOkm are common, as shown in Fig. 7. 
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Surfdc? Mean Temperature Structure For Seotember 
1S77 (From Renner, 1977). 
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IV 



EXPEaiMENT 



A. CEJECTIVE 

For ocean monitoring from space to become effective, the 
interior metiers of the oceans must be inferred from the 
resulting signature of surface properties. The ability to 
define the subsurface structure accurately by this method is 
becoming increasingly utilized and understood. 

legeckis and Gordon (1982) described the ability of sat- 
ellite infrarsd imagery tc locate sea surface temperature 
fronts, and the relation cf these surface features to sub- 
surface tempaiaturs structure. For the Brazil Current and 
two warm cere eddies, they noted the mixed layer depth cor- 
related with the surface temperature patterns. 

Bernstein, Breaker and Whrit ner (1977) examined the eddy 
formations in the California Current from combined ship, air 
and satellite platforms and found that the sea surface 
temperature patterns conformed to subsurface warm and cold 
mesoscale variations. They concluded that in unstable 
regions such as the California Current, surface and main 
thermccline temperature distributions have relation ships to 
each other. Thus, in such a region, remotely sensed sea 
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surface temperature patterrs would be able zo interpret sub- 
surface circulation. 

From these recent examples of sareilire-determined sub- 
surface structure and from other oceanographic literature, 
it appears that the regicns of the world's oceans where 
large instabilities occur (i. e. , boundary currents, coastal 
upwelling zones, etc.), the subsurface dynamic features are 
readily deduced from the surface observations. Areas of the 
seas which are quiescent and usually located in the central 
oceans have net been as easy to describe accurately. In 
these regicns, the use cf either dynamical or empirical 
methods cculd be the only means of forecasting the thermal 
structure cf the oceans. 

Numerical models are based on dynamical principles and 
will ultimately provide the basis for predicting oceanic 
processes. Empirical methods have been determined in the 
past for certain regions and for limited time periods, but 
do not appear tc be valid universally. If relationships 
could be devised for certain water masses or oceanic 
regimes, then the possibility exists of piecing together the 
subsurface thermal features in the different regions by uti- 
lizing remote sensing to define certain surface 

measurements . 
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Toward rhis goal, the first objective of this thesis is 
zo determine whar particular relationships exist between sea 
surface temperature and subsurface thermal structure for a 
relatively "quiet" region in the northeast Pacific Ocean. 
By examining the climatology, an empirical set of equations 
could then be derived between sea surface temperature, mixed 
layer depth, thermocline gradient and the depths of promi- 
nent points within the structure. 

The second objective is to produce a predicted vertical 
temperature profile from the best fit regression equations 
for a particular water region and SST at a future time and 
different location. When these profiles are compared to 
actual temperature- de pth profiles, accuracy of this statis- 
tical method will be determined. 

B. CfiTA ACQUISITION 

Curing the period 5 to 14 September 1977, the U.S. Naval 
oceanographic survey vessel USNS SILAS BENT conducted ocean- 
ographic surveys to determine the thsrmohaline and sound 
velocity structures across the Califiornia Current, the 
Transition Zone and Subtropical Water. For this thesis, tne 
only data utilized were along a 1155 km meridional track 
south of Ocean Weather Station "P", as displayed in Fig. 8. 
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The data ccnsisted of 28 vertical temperature profiles 
spaced appro ximaxely 42 km apart; they spanned the period 9 
to 11 September 1977. 

Weather conditions during this period were governed by a 
weak staticnary low pressure system with mostly low clcud 
cover and periods of light drizzle and showers. The wind 
speed varied from 5 to 20 knots from several directions 
causing wave heights of approximately 3 to 6m along the 
en+ire track (FNCC, 1977). 

Expendable bathythermographs (XBT*s) were used to define 
the thermal structure. The resulting temperature traces 
plotted onboard ship were coded on bathythermograph log 
sheets, which were utilized as the data points for recon- 
struction cf the thermal profiles. Characteristic of this 
procedure is the loss of the fine thermal structure, as only 
prcminent points cn the temperature profile are recorded to 
give the general shape and location of important features. 

Another problem associated with this method of data 
acquisition is the varying probability that the observer on 
the vessel at the time of the X ET drop correctly chose the 
proper points to reflect adequately the mixed layer depth, 
thermocline gradient and temperature inversions. Incorrect 
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Figure 8. Lccaricn of XBl DATA Curing 3-11 Sepreabej 



1977. 
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recording cf the values along with maif uncrioning in 
10% cf the XEls can also te a source of potential error. 

C. ASSESSKENT CF TEMFERAIUEE STRUCTORE 

Tc begin the analysis cf subsurface thermal structure, 
realizing the “crude" nature of the data, the temperature- 
depth points were plotted cn a computer graphics system and 
the 28 thermal profiles were reconstructed. From these pre- 
files a vertical cross-section of the meridional temperature 
structure was drawn. Fig. S. It was noted immediately that 
numerous temperature inversions existed along the northern 
and central regiens of the track, but they ceased at approx- 
imately 41N. This subsurface border clearly defines the 
Subarctic Front which divides the Subarctic from Subtropic 
Maters; while at the surface there is no indication cf the 
division. The temperature inversions varied in strength, 
but generally shewed an increase in depth in a southward 
direction . 

Ihe upper mixed layer displayed little variablility with 
few transients found within the structure. There were no 
large eddies rcted in the surface layer and the mixed layer 
depth appeared relatively constant along the entire track. 
Weak mescscale features are probably present, but the 
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Figure 9. Meridicnal Cr oss-sect ion of Temperature at 145W. 
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spatial r-sscluxicn prohibits their detection. The thermo- 
cline was also fairly constant in depth with the strongest 
gradients between 30 to 6C m. The isotherms in the lower 
layers characteristically deepened to the south, with sev- 
eral areas of large variability located near the temperature 
inversions. The large warm intrusion located at U7N may be 
caused by errors in the data, as a structure this size in 
the Subarctic Region does not appear to be congruent with 
the ether features, and is not mentioned in the research 
li "sratur 9 . 

To assess further the vertical thermal structure, promi- 
nent points, gradients and zones were defined on a typical 
temperature profile from the test region. Fig. 10 illus- 
trates the variables that were used to analyze sta-ist ically 
the thermal profile. Descriptions of these variables are: 



• SST - Sea surface temperature was defined as the sur- 
face temperature reported on the bat hyt her moqraoh logs. 
Since this is the temperature recorded on the XBt trace 
when the XBT orobe entered the water, a slight dispar- 
ity from remotely-sensed SST may occur due to the 
response time of the XBT system. 



• dLD - Mixed layer depth was visually determined as the 
greatest depth where the isothermal or slightly 
negative temperature gradient in the surface layer 
changed to the stronger gradient of the underlying 
ther mccline . 






BCT - Bottom of seasonal thermocline was 
determined as the depth where the seasonal t 
markedly changes into a weaker negative t 
gradient below. 



visually 
her mo cl in 
emperat ur 
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• TBOI - Temperature at the bottom of the seasonal ther- 
mocline. 

• DT - Difference in temperature between the temperature 
at the mixed layer deptn and at the bottom of seasonal 
ther mccline . 

• DZ - Difference in depth between the mixed layer depth 
and bottom of seasonal ther mocline. 

• EOZ - Depth of the bottom of the main thermocline was 
visually determined as the depth where the thermocline 
changes to an isothermal or positive temperature 
gradient *ust above the temperature inversion region. 

• TBUZ - Temperature at the bottom of the main thermo- 
cline . 

• DTBU2 - Difference in temperature between the 
temperature at the mixed layer depth and at the bottom 
of main thermocline. 

• DZBUZ - Difference in deoth between the mixed layer 
depth and the bottom of main thermocline. 

• T460 - Temperature at reference deptn of 460 m which 
was selected arbitrarily at maximum depth of XBT 
reports. 

• DT5 - Temperature change in upper 5m of thermocline 
below the mixed layer depth. 

• DT15 - Temperature channge in upper 15 m of thermocline 
below the mixed layer depth. 

• DT50 - Temperature change in upper 50 m of thermocline 
below the mixed layer depth. 



DT100 - Temperature chance in upper 1 00 m of ther 
dine below has mixed layer depth. 

DIS - Distance from northern most XBT drop. 



mc- 
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Figure 10. Model of Vertical Temperature Structure. 
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This assessment defines the general structure primarily 
in the upper layers above the permanent halocline. The 
detail of the smaller structure such as transients in the 
mixed layer and multiple thermoclines are not included 
because they are variable and minor with respect to the 
principal structure. The temperature inversion feature 
below the main thermocline was not charact erir ed as its 
occurrence in the region varies in location and intensity 
during the test period. In those zones within the Subarctic 
Segicn where a sub-thermccline temperature maximum exists 
permanently during the heating season, additional variables 
should be used to define this part of the thermal profile. 

Initially, the sea surface temperature gradient was used 
as a variable in the analysis. However, problems occurred 
in computing consistent gradients along the trade as the 
calculated values were unrealistic and inaccurate compared 
with actual conditions at the time; this variable was 
dropped from this study. 

These verities provided the statistical base from which 
all ccrrelaticns and regression equations were derived. The 
values could then distinguish the fluctuating and transient 
features from the median thermal structure. Tables I, II 
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and III list tha values from the analyzed thermal profiles 



by XET drop number. 



TABLE I 

locations and Times of XBT Measurements 



XBT 


CATE^IME(Z) 


LAT (N) 


LONG ( W) 


DISTANCE (km) 


339 


0906 13 


SEP 


77 


47. 98 


144. 98 


0 


240 


0SC804 


SEP 


77 


47. 65 


144. 93 


36 


341 


091003 


SEP 


77 


47. 20 


144 . 98 


86 


342 


CS1200 


SEP 


77 


46. 82 


144. 32 


128 


343 


091800 


SEP 


77 


46. 07 


145. 05 


212 


244 


0S2000 


SEP 


77 


45. 70 


145. 04 


253 


345 


092205 


SEP 


77 


45. 33 


145.04 


294 


346 


100000 


SEP 


77 


45. 02 


145.04 


329 


347 


100308 


SEP 


77 


44. 73 


145. 00 


361 


248 


1005C8 


SEP 


77 


44. 37 


144. 98 


401 


349 


100708 


SEP 


77 


44. 00 


144. 98 


442 


250 


1CC904 


SEP 


77 


43. 63 


144. 98 


483 


351 


101100 


SEP 


77 


43. 28 


144 . 97 


522 


252 


1C1400 


SEP 


77 


43. 07 


144 . 95 


546 


353 


101600 


SEP 


77 


42. 70 


144. 93 


587 


354 


1C1800 


SEP 


77 


42. 32 


144 . 95 


629 


355 


102000 


SEP 


77 


41. 93 


144.93 


672 


256 


110002 


SEP 


77 


41. 55 


145. 00 


714 


357 


110203 


SEP 


77 


41. 17 


145. 00 


757 


258 


1 104C6 


SEP 


77 


40. 75 


145. 00 


803 


359 


1 10604 


SEP 


77 


40. 37 


145 . 00 


846 


360 


1 1 1000 


SEP 


77 


39. 97 


144. 98 


890 


361 


1 1 1200 


SEP 


77 


39. 60 


145. 10 


931 


362 


1 1 1400 


SEP 


77 


39. 22 


145. 25 


973 


363 


1 11600 


SEP 


77 


38. 83 


145.37 


1017 


264 


1 12304 


SEP 


77 


36. 30 


145.47 


1076 


365 


120107 


SEP 


77 


37. 93 


145. 53 


1117 


266 


1203C7 


SEP 


77 


37. 58 


145.67 


1156 
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XBT 

339 

340 

341 

342 

343 

344 

3 45 

3 46 

3 47 

348 

349 

350 

351 

352 

353 

354 

355 

3 56 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 



TABLE II 

Values of SET, MLD, Thermoc lir. e Gradients and T460 
SST (C) t!ID (m) DT5(c) DT15(c) DT50(C) DTIOO(c) T460 (c) 



13.4 




1.5 


4.6 


13 . 6 


-24 


0.3 


1,5 


14. 1 


- 25 


2.5 


4.3 


14.4 


-25 


2, 3 


5,6 


15. 1 


-19 


2.4 


4.8 


15.3 


-32 


0. 3 


4.6 


15.6 


-29 


3. 1 


4,9 


15.7 


-19 


0. 4 


4.7 


16,6 


-27 


0.2 


4.3 


16.2 


-20 


0. 5 


4.7 


17.0 


-21 


0.4 


5,3 


18.1 


-20 


0, 2 


6,2 


18. 3 


-30 


3,6 


4.2 


18.1 


-20 


0.3 


5,8 


18.4 


-20 


0. 7 


4.0 


18.1 


-20 


0. 6 


5.4 


19.3 


-21 


0.6 


5.4 


19.7 


-15 


0. 2 


6.9 


19.5 


-23 


0.6 


6.4 


20. 1 


-16 


0. 3 


2. 1 


20.3 


-22 


0. 8 


7.0 


20. 1 


-16 


0. 5 


2.6 


20.5 


-30 


2.4 


6. 1 


21.5 


-22 


0. 3 


5.7 


21.4 


-32 


3. 8 


6.4 


22.5 


-27 


0. 4 


6. 2 


22.3 


-22 


0,4 


5.6 


22.7 


-20 


3. 3 


5.2 



7.1 


6.6 


4. 0 


6.6 


7.4 


4.1 


6.0 


6. 2 


5.0 


7.1 


7.6 


4.2 


7.3 


8. 1 


4.6 


7.7 


7.4 


4 .9 


7.3 


7. 9 


6.2 


6.7 


7. 9 


4.7 


7.9 


9. 1 


4.7 


7.4 


8.4 


4.8 


8.0 


8. 5 


4.9 


9.1 


10.3 


4.9 


7.8 


10. 2 


5.0 


8.4 


8.7 


5. 1 


9.0 


9.9 


5.0 


8.6 


9.3 


5.2 


9.3 


10.3 


5.3 


9.6 


10.8 


5.3 


9.7 


10. 7 


5.3 


9.4 


10.3 


5.6 


10.0 


10.8 


5.6 


8.3 


10.3 


5.9 


9.3 


10.4 


5 . 5 


10.5 


11.6 


5.8 


9.8 


11.4 


5.8 


9.8 


11.1 


6.5 


7.6 


9.6 


6.9 


8. 1 


10.4 


6.8 
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TABLE III 



Values of Thermocline Temperatures and Deorhs 



111 


EOT (tn) 


IIO iQ) 


R1 (c) 


12 (m) 


BUZ (m) 


TBUZ (c) 


DTBUZfc) 


DZBUZ 


339 


-88 


5.4 


7.4 


55 


-95 


5. 3 


7. 2 


62 


340 


-30 


6.5 


6.8 


56 


-109 


5.9 


7.4 


85 


34 1 


-70 


7.6 


5.9 


45 


-110 


7.4 


6. 3 


85 


34 2 


-79 


6.6 


7.2 


54 


-128 


6.2 


7.6 


103 


34 3 


-81 


6.6 


8.0 


62 


-138 


6 . 5 


8. 1 


1 19 


344 


-81 


7. 1 


7.7 


49 


-112 


6.7 


8. 1 


80 


345 


-68 


7.8 


7.1 


39 


-1 1 2 


7.0 


7. 9 


83 


34 6 


-57 


8.4 


6.8 


38 


-100 


7. 3 


7.9 


81 


34 7 


-68 


8.6 


7.8 


41 


-114 


7. 3 


9. 1 


87 


348 


-65 


8.4 


7.3 


45 


-143 


7.6 


8. 1 


128 


34 9 


-50 


9.1 


7.5 


29 


-109 


8. 1 


8. 5 


88 


350 


-56 


8.8 


9 .0 


36 


-123 


7.5 


10.3 


103 


35 1 


-57 


10.7 


7.0 


27 


-148 


7. 2 


10. 5 


1 18 


35 2 


-73 


9.0 


8.5 


53 


-158 


7. 3 


9.7 


138 


353 


-67 


9.4 


8.8 


47 


-132 


7.7 


10. 5 


112 


354 


-63 


9.3 


8.6 


43 


-120 


8.6 


9.3 


100 


35 5 


-59 


9.9 


9.1 


38 


-130 


8.6 


10. 4 


109 


356 


-54 


9.9 


9.5 


39 


-134 


8.5 


10. 9 


1 19 


35 7 


-62 


10.1 


9.1 


39 


-127 


8.4 


10. 8 


104 


35 8 


-50 


11.3 


8.6 


34 


-120 


9- 3 


10.6 


104 


35 9 


-53 


10.8 


9.4 


31 


-95 


9.4 
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To examine the variations and trends of the individual 
therial profile features along the track, graphs of the 
variables as related to distance and latitude are plotted in 
Appendix A. These diagrams illustrate some of the concepts 
that were discussed in previous chapters pertaining to the 
variability in thermal structure within a particular region. 

The sea surface temperature over the entire track as 
shewn in Fig. A.1 increases linearly to the south where the 
mean gradient was 0.8C/100km. The strongest gradient 
occurred in the vicinity of h 4N and reached a value of 
2. 3C/100km. Fig. 11 shews the sea surface temperature in 
the North Pacific as compiled from ship reports and clima- 
tology by Fleet Numerical Cceanography Center for 10 Septem- 
ber 1977. In the region of data acquisition, the meridional 
temperature gradient and surface isotherm structure closely 
compare tc the surface temperature versus distance plot 
(Fig. A.1). Satellite infrared imagery during this time 
period was obstructed due tc the thick cloud cover over the 
area; therefere, this method of sea surface temperature 
determination was not used as was desired (Tabata and 
Kimber, 1 979) . 
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Figure 11. Pacific Ocean Sea Surface Temperature on 10 
September 1977 (Fleer Numerical Oceanography 
Center, 1 977) . 
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Fig. A. 2 shews ehe dspth of the mixed layer along the 
track. The mean dspth was -22. 9 ra with a maximum and mini- 
mum depth of -32 and -15 m, respectively. The layer's gen- 
eral trend of a slight shallowing by about 5 m in the 
southward direction is observed until 40N is reached where 
the trend reverses and becomes deeper. This location of 
change in the mixed layer depth trend accurately coincides 
with the Subarctic Front. 

The four thermocline gradient variables in Figs. A. 3 
through A. 6 displayed a measurable trend along the track. 
The variability of the temperature change in the upper 5 m 
of the thermocline (DT5) was large in the northern and 
southern ends cf the track. In the center, with the excep- 
tion of cne data point, a mean value of 0.3C/5m was 
observed. The boundaries of these three zones of the BT5 
variable agreed with the division of the track into the 
three regions, i.e. , Subarctic Water, Transitional Zone and 
Subtropic Water. This effect was noted in many of the ther- 
mal structure variable plots. DT15, DT50 and DT100 all 
increased with decreasing latitude, but the general trends 
differed between each variable. DT15 showed the largest 
range of variability with only a slight increase in 
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temperature gradient, while DT50 and DT100 borh increased at 
a greater rate. When plcrted collectively, as shown in ?ig. 
12, xhe relationships become apparenr. Generally, with 
increasing sea surface temperature, the negative temperature 
gradient of the thermccline increases slowly at a relatively 
constant rate, 

Iwc ether depth-dependent variables on the vertical 
temperature profile, besides the mixed layer depth, are the 
bottom of the seasonal and main thermoclinas , BOT and BUZ, 
respectively. Figs, A, 7 and A, 1 1 display these variables as 
they vary with latitude. The variability of the plots may 
be due in part tc internal waves with typical amplitudes of 
15 m for the ECT and 30 m for the BUZ, 

The BOT’s general trend southward along the track to 
about 40N shewed a slight decrease in depth from -80 to 
-60 m. Again, this location coincides with the Subarctic 
Front as did the MLD, The BOT then increases sharply as the 
trend reverses tc a depth of -90 tn. 

For the BUZ, the trend is just the opposite of the EOT 
in the Subarctic Water and Transitional Zone, The depth of 
the bottom of the main thermocline increases slowly until 
reaching the Subarctic Front, where the maximum variation in 
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Ccmparison of DT5 , DT15, DT50 and DT100 Versus 
Latitude. 
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dapth is apprcximat aly 85 m. In Subtropic Warer, the rrsnd 
of the BDZ increases in depth at a greater rate along the 
southern and cf the track. 

When the KID, BOT and EUZ are plotted together, as shown 
in Fig. 13, additional relationships appear. In the Subarc- 
tic and Transition Regions, north of 40N, the thickness cf 
the upper layers above the bottom of the seasonal thermo- 
cline (BOT) decreased in a southward direction. With the 
MLD trend at a relatively constant depth, the only thermal 
layer to change in magnitude appreciably is the thermocline. 
This suggests that the seasonal thermocline, with its 
negative temperature gradient, becomes stronger as the SST 
increases and as the thermocline thickness decreases. 
Indeed, the temperature difference (DT) and thickness (DZ) 
of the seasonal thermocline * displayed in Figs. A. 9 and 
A. 10, respectively, agree with this suggestion. 

As previously mentioned, the dep+h of the bottom cf the 
main thermocline gradually deepens toward the south along 
the track. Figs. A. 13 and A. 14 show the increasing differ- 
ence in temperature and depth from the mixed layer depth, 
DTBUZ and DZBOZ, respectively. These tendencies, in con- 

junction with the trend of the BOT, indicate that the 
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negativa temperature gradient at the lower part of the main 
tharmocline becomes weaker and thicker reward the south in 
the Subarctic and Transiticnal Regions. 

For Subtropical Water, south of 40N, the BOT has 

reversed its trend at the Subarctic Front and increases in 

depth. The EDZ continues to deepen, but at a greater rats 

than to the north. Temperatures at these two locations on 
the thermal profile, I30T and TBUZ, slowly increase linearly 
with decreasing latitude at nearly the same rate of 
0.5C/100km. These values suggest that the negative 
temperature gradient at the lower part of the main thermo- 
cline continues to weaken and thicken, but at a deeper depth 
than in the northern waters. At the upper layers, the 

parameters of the seasonal thermocline, DT and DZ, increase 
with increasing SST. 

D. STATISTIC?! EBOCECURE 

To predict the subsurface thermal structure from a given 
sea surface temperature, a statistical approach was utilized 
to determine the relationships between the thermal structure 
variables. The data analysis was done using the BIOMED sta- 
tistical computer programs on the Naval Postgraduate 
School’s lEM 3033 computer system. Details of the BIOMED 
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computer pregrams used in the analysis are provided in 
Appendix E. Basic statistics (mean and standard devia- 
tions), scatter diagrams, correlation and linear regression- 
analysis were performed on all variables. 

As part cf the regression analysis, the correlation 
coefficients fer all combinations of pairs of the sixteen 
variables were computed. These simple correlations are an 
estimate cf the strength cf the linear association between 
two variables and do not imply the cause for the relaticn- 
ship cr which variables affect the variation of the other. 

Scatter diagrams were used in conjunction with the 
regression analysis to assist in refining the correlations 
between variables. By visually removing poor data points 
that detracted from trends within specific regions, the 
analysis cculd be improved. After determining the best fit 
correlations, least square linear regression equations are 
computed between variables and also with sea surface 
temperature as the independent variable. This enabled a 
vertical temperature profile to be determined by calculating 
the values fer the variables with the only input being the 
sea surface temperature. 



62 



The values cf the 16 thermal structure variables from 
the 28 temperature profiles along the track were correlated 
for all ccmbications with the results shown in Table IV. 
Several strong relationships appear with the correlation 
coefficients cf a few variables being greater than 0.750. 

The sea surface temperature had the greatest number of 
strong correlations which greatly enhances the idea of pre- 
dicticn cf the thermal prcfile from the SST. The highest 
correlaticn coefficient (0.993) occurred when the SST was 
compared with the distance along the track. This was 
expected due tc the relatively constant north-south surface 
temperature gradient in this region of the Northeast 
Pa ci f ic. 

ether strong correlations with the sea surface 
temperature were the temperatures at the bottom of the sea- 
sonal and main thermocline, T30T (0.944) and T3UZ (0.949), 
respectively. When plotted together, as shown in Fig. 14, 
the trends clearly appear. The variables' similar relation- 
ships imply that the entire vertical temperature profile 
"leans" reward warmer temperatures as the small negative 
thermal gradient in the deeper layers increase with increas- 
ing sea surface temperature. The change in temperature in 
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Corrsiaticn Coafficients of Thermal Variables 
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the upper 50 and 100 m cf the thermocline also displayed 
strong direct relationships with the SST. 

Two pairs of variables which exhibited strong negative 
correlaticns , and for similar reasons, were BOT-DZ (-0.949) 
and EUZ-DZEOZ (-0.983). These combinations indicate that a 
relationship exists between both the BOT and BUZ and the 
MLD. This is due to the relatively small variations in the 
mixed layer depth when compared to the large fluctuations of 
the EOT and BGZ. Therefore, as the bottom cf the seasonal 
and main thermoclines varied, the depth difference between 
these variables and the MLD similarly changed. This could 
provide an important link in determining the mixed layer 
depth if a viable connection cculd be established between 
the EOT or BUZ and SST. However, according to the correla- 
tion coefficients comparing the SST to MLD (0.224), to BOT 
(-0. 203) and to BUZ (-0.682), the strength of these rela- 
tionships was weak as shown in Pigs. 15, 16, and 17. Corre- 
laticns of the ether variables with the mixed layer depth 
revealed little as to possible methods to determine this 
important thermal structure variable. 

Due to the noticeable differences in thermal structure 
characteristics among the various regions encountered alcng 
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Latitude . 



66 




Figure 15. Cc iipc.rison of MLD and SST Versus Latitude 
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the track, the data ware next divided into Subarctic, Tran- 
sitional, Frontal and Subtropical groupings. These divi- 
sions were detsririned froir the previously discussed regional 
characteristics and thermal analysis of Chapters II and III, 
and in retrospect, they proved to be an accurate reflection 
of the boundary locations at the time of data acquisition. 

The Subarctic Water characteristics were traced from the 
northern-most vertical temperature profile on the track tc 
XBT station number 345 at 45.33N. The Transitional Zone was 
observed from XET station 346 to station 357 at 41.17N fol- 
lowed by the Subarctic Front between stations 358 and 359. 
Beginning at XBT station 360 at (39.97N), and continuing 
south to the end of the track, was an ocean region contain- 
ing properties of the Subtropic Hater. Fig. 18 illustrates 
the divisions of the track data points into the various 
regions. 

To observe the separation of the data by groupings bet- 
ter, and to refine further the correlation coe f f icients , 
scatter diagrams of several of the thermal variables ware 
plotted using the BIOHED computer programs. Since the only 
variables that correlated strongly with SST were other 
temperature-dependent variables such as TBOT, DT, TBUZ and 
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DTBOZ, a closer inspecticn of the depth- d spend enr variables 
was desired. To ascertain whether the sea surface 

temperarure can define the depth of cerrain variables is 
necessary in the development of the nhermal profile. 
Therefore, regional scatxer diagrams for the xhrae primary 
depth dependent variables (MLD, BOT, BUZ) and SST as com- 
pared zo the ether variables were plotxed and analyzed. 

The scatxer diagrams for sea surface temperaxure when 
compared to the ether variables are lisxed in Appendix C. 
Due to the high correlaticn coefficients of these thermal 
variables the trends of the data points were clearly visi- 
ble. Obvious linear relations cccur between SST and DT50, 
DT100, DT, TECT, TBUZ and DTBDZ through the various regions. 
DT5 and DT15 also displayed an element of linearity when 
several of the mere distant data points ware removed. The 
computer-derived correlation coefficients for these two 
variables were excessively low due to the remote points. 
For the variable DT5 shown in Fig. C. 2, the data behaved in 
a strong linear fashion in the Transitional Region as a 
temperature change of 0.3C in the upper 5 m of the xhermc- 
cline was the mean. All cf the thermoclina gradient vari- 
ables, DT5, DT15, DT50 and DT100, could then be modeled by a 
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linear eguaticn related tc the SST, at least in the Tran- 
sition Region. 

Two scatter diagrams which compare SST to BUZ and DZBUZ, 
Figs. C.IC and C.13, illustrate a perceivable linear trend 
between the variables while the correlation coefficients 
reflect a moderate level of association, (-0.687) and 
(0.678), respectively. Again, with several of the remote 
data points removed, an improved correlation coefficient 
could be computed which better portrays the relationship 
with SST. 

ether depth-dependent variables, for which a visual 
relationship was observed with SST in the scatter diagrams, 
were for ECT and DZ, Figs. C.6 and C.9. Both exhibited a 
"V* or wedge shape with data for the Subarctic and Subtropi- 
cal Haters on either wing and data for the Transitional zone 
in the center. If the northern or southern regional data 
were removed, then a higher correlation could be obtained 
utilizing that from the Transitional Region alone. However, 
with the visual trend in this region closely following the 
trend of data from the Subarctic Hater, data from the north- 
ern end of the track could be used to assist in the accurate 
determination of relationships with the SST. The 
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CO mbinaricn of data from thass two regions is raasonahla 
since the Iccaticr. of the Subarctic Front is on the southern 
boundary cf the Transition Zone, and the waner characteris- 
tics are siirilax. Therefore, the Subrropic Water data 
points should be removed tc derive the correlations with 
data from only the two northern regions. 

The twc data points in the regional scatter diagrams for 
the Subarctic Front consistently followed the Transitional 
Zone data trends. As shewn in Appendix C, for the thermal 
variables EOT, DZ, BUZ and DZEUZ, the Front provided the 
sharp break between the Transitional and Subtropical 
Regions, but in the scatter plots for the remaining thermal 
variables the Subarctic Front fell within the trend of the 
twc northern regions. 

The scatter plot comparing mixed layer depth tc SSI, 
shewn in Fig. C.l, initially appears not to reveal any pat- 
tern which cculd be used to infer a relationship with the 
SST. Hewever, if the Subtropical Water data points are not 
included, a trend of decreasing depth is noticed and the 
correlation ccefficient value of 0.224 between these two 
variables should increase appreciably. 
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The regicnal scatter diagrams comparing mixed layer 
depth to the ether variables did not exhibit any s-^rong cor- 
relations and were not included. However, several weak pat- 
terns were observed which indicated a trend of the MLD 
slowly decreasing in depth in a southward direction for the 
northern twe regions. All scatter plots displayed large 
variability among the data points which agrees with the lew 
correlation coefficients for the variable. 

Appendices D and E contain the regional scatter diagrams 
for ECT and EUZ, respectively. These plots show a greater 
linear dependence among the variables, and therefore higher 
correlation coefficients, than was shown by the mixed layer 
depth. In most cases when the Subtropical Water data are 
removed, leaving only data from the Subarctic, Transition 
and Frontal regions, a stronger relationship was observed. 
Two noteworthy scatter diagrams which had strong 
negative correlations were BOT-DZ and BUZ-DZBUZ as shown in 
Figs. 19 and 20. As previously noted, these two variables 
with high correlation coefficients could provide the neces- 
sary link to determine the mixed layer depth. Since DZ and 
DZBUZ are the differences in depth from the HLD to the EOT 
and EUZ, respectively, a pair of simple equations define the 
relationship. 
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MLD 



EOT ♦ DZ 



MLD = EDZ + CZBUZ 

note: MLC, EOT and BUZ are negative values (depth) 

The only weak "link" in the two equations is the deter- 
mination of the EOT and BGZ directly from SST. To improve 
the correlation coefficients for these two variables as sug- 
gested by the regional scatter diagrams, data from the Sub- 
tropical Hater and several distant data points which 
detracted from the established trends were removed. In sev- 
eral cf the plots, primarily fcr the temperature-dependent 
variables, where the trend was actually stronger and more 
accurate with data from the Subtropical Region included, the 
correlation ooafficients were improved. The correlation 
coefficients fcr the best-fit trends are listed in Table V. 
Indications are made to show which coefficients increased or 
decreased from the values presented in Table IV. 

An example cf a correlation coefficient that decreased 
due to the removal of data points and yet provided a 
stronger trend was the KID versus DT5 combination which 
changed from -0.484 to -0.197. In this case, the majority 



76 



DEPTH OF BOTTOM OF SEASONAL THERMOCUNE <M) 
-no -100 -90 -«0 -70 -60 -90 




A 



X 




A 



LEGEND 
• ■ SUSAECriC 
A - TRAi’^LSmOMAL 
X- FRONT 
a - SUBTROPICAL 






□ 






□ 




DZ - DEPTH DIFFERENCE FROM MLD TO BOT <M) 



Figur-s 19. 



Recicna 1 



Scatter Diagram 



DZ V'^rsus BOT. 



77 



DEPTH OP BOTTOM OF MAIN THERMOCUNE CM) 

-a»o -1*0 -170 -ISO -»Jo -no -*o 



I 



• • 






L X 
A 

a ^ 



4 



LEGEND 
• - SUSASCTIC 
A - TRANSITIONAL 
X - FRONT 
am SUBTROPICAL 



Bm 

00 



70 *0 no 130 190 

DZSUZ - DEPTH DIFFERENCE FBOM MLD TO BUZ CMJ 



170 



1*0 



Figure 20. 



Regicral Scat 



er DiagraoQ 



DZBUZ Versus BUZ. 



78 



TABLE V 



B<=st-Fit Ccrr€laticn Ccefficients of Tharmal Variables 
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of data points were grouped about the DT5 value of 0.4C per 
5 m and a MLE value of -23 m . The attempt to capture the 
trend of the Transition Region was made by excluding the 
distant scattered points; however, this detracted from the 
overall trend and reduced the correlation coefficient, 

A large number of correlation coefficients improved sig- 
nificantly by this method of data exclusion. Most notice- 
able is the variable SST, which, (except for the combination 
with DT5) increased all other coefficients to greater than 
0.645 with six increasing above 0.900. The important corre- 
lations of SST with EOT and BUZ improved from -0,203 and 
-0.682 to 0.776 and -0.790, respectively. The correlation 
coefficient of SST versus EOT even reversed sign. 

This increased strength in the various correlations 
between variables provided the confidence to infer that usa- 
ble relationships exist in the vertical thermal profile. 
Since the variable SST had the greatest number of strong 
correlations with the other variables, the statistical 
method of deriving linear regression equations using only 
the input of SST appeared practical. 
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E. EEGRESSICN ECUATICN ANALYSIS 

From the kest fit correlation coefficients of the vari 
able SST in ccmbination with the other variables, the fcl 
lowing linear regression equations were derived: 

MLD = (1.2C01) SSI - 43.300 

CT5 = (0.03629) SST - 0.21573 

DT15 = (0. 16058) SST + 2.2271 

DT50 = (0.48462) SST - 0.14105 

DT100 = (0.48506) SST + 0. 52999 

BCT = (3.9611) SST - 133.08 

TBOT = (0.61762) SST - 1.9317 

DT = (0.41 729) SST + 0.94677 

DZ = (-2.7611) SST + 89.782 

BUZ = (-7.5247) SST - 0.00793 

TEUZ = (0.4851 0) SST - 0.66001 

DTEUZ = (0.55345) SST - 0.42130 

DZEUZ = (7.7697) SST -27.820 
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To construct a complete vertical temperature profile, 
since the variables previously mentioned only defined the 
upper domain cf the theriral structure, a new variable was 
needed from Icwer depths. The temperature at 460 m was cho- 
sen because it was the lowest datum point available from the 
XBT data provided along the track. This type of variable 
was net particularly desired as it relates a temperature to 
a specific depth rather than to a feature on the thermal 
profile. A better variable would be at a reference depth 
where the temperature gradient becomes isothermal in the 
deeper regions of *he ocean. However, this was not availa- 
ble from the limited depth of the data and the T460 variable 
was established. The same methods of analysis were employed 
on this variable as were conducted on the others. When cem- 
oarsd to the sea surface temperature, a correlation coeffi- 
cient of 0.331 was computed along with a strong linear 
trend in the regional scatter diagram. Therefore, a linear 
regression equation could be used to determine the bottom 
point of the predicted temperature profiles at 460 me~ers 
from the variable SST. The linear regression equation for 
T460 was computed as: 

T460 = (0.21334) SST + 1.3109 
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Two ether iaiccrtant relationships which were derived by 
regression analysis, but they do net use SST as the indepen- 
dent variable: 



DZ = (-0.85394) BCT - 12.693 

DZBUZ = (-1.0359) BUZ - 27.392 

These equations use the computed values of BOT and BUZ 
as determined from the SST to calculate the thickness cf the 
seasonal and main thermoclines . The magnitude of these 
thicknesses (CZ, DZEUZ) when added to the BOT or BUZ, 
respectively, can then define the mixed layer depth. 

Due to the relatively low correlation coefficient 

between SST and MLD for the SST-derived MLD equation 
(0.535), possibly a combination of equations which define 
3LD with higher correlations could improve the accuracy in 
determining the mixed layer depth. The regression equation 
for the SST-derived BUZ, with a correlation coefficient of 
-0.790 between the two variables, along with the BDZ-derived 
DZBUZ manticned above with a correlation coefficient of 
-0.983, seemingly provide the strongest combination. 
Therefore, the group of equations used for the mixed layer 
depth determination were: 
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BUZ 



(-7.5247) SST - 0.30793 



DZBOZ = (- 1.0359) BUZ - 27.892 
MLD = BUZ + DZBUZ 

Ihres ether combinatiens of regression eguarions could 
have been used to determine the mixed layer depth# but were 
not chosen because of their lower correlarion coefficients 
between variables. The equations and correlation coeffi- 
cients (r) are listed belcw: 



BOT = (3.S611) 


SST - 


133.08 


r = 0.7 76 


DZ = (-2.7611) 


SST + 


89. 732 


r = -0.645 


MLD = BCT + CZ 









EOT = 


(3.S6 11) SST - 


133.08 


r = 0.776 


DZ = 


(-0.65394) EOT 


- 12.693 


r = -0.94 



MLD = BCT ♦ DZ 
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EQZ 



(-7.5247) SST - 0. 007 93 



-0. 790 



DZBUZ 



(7.7697) SST - 27.820 



MLD = BUZ + DZBUZ 



r = 0.766 



With ths aiD determined, the “empararure at this depth 
was needed tc provide the fundamental tamparature from which 
the tharmccline could be computed. A variable TMLD was 
added to represent the temparatura at the mixed layer depth 
and was defined as the sum of TEUZ and DTBUZ. Since both 
tamperature-dependent variables had very high correlations 
with the SST, 0.949 and 0.960 respectively, strong confi- 
dence could be placed in the resulting variable, TMLD. 

An additional feature in calculating TMLD by this method 
was to allow for a slight negative temperature gradient tc 
occur in the mixed layer. Many ocean thermal structure mod- 
els assume that this upper layer will be isothermal when 
actually a realistic profile often includes a temperature 
gradient . 

The four thermal gradient variables DT5, DT15, DT50 and 
DT10C, as determined from the SST, were subtracted from the 
temperature and depth of the mixed layer to form the thermo- 
cline. This structure, along with the SSI and mixed layer 
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depth values, provided the upper portion of the ver~ical 
temperature profile. 

In the region below the seasonal rhermccline the 
temperature-dependent variables, TEOT and IBUZ, and depth- 
dependent variables, EOT and BDZ, were computed directly 
from the input of SST and applied to the temperature pro- 
file. Several times during the plotting of these points in 
the testing and ccmpariscr phase of the thesis, the Icca-^ion 
of the DT100 variable would occur below both the BOX and BUZ 
or between them. In these cases it was simple to adjust the 
sequence of points to allow the proper construction of the 
temperature profile. 

Finally, in the generation of the temperature profiles, 
the variable T46C was used to provide the lower base of the 
structure. Since T460 was an artificially-generated vari- 
able, defined at a constant depth, the characteristic small 
temperature variations at that depth should still be accu- 
rately reflected. 

A total of nine points were used to define the entire 
vertical temperature profile. These were plotted on a 
TSKT5CNIX 618 graphic display system and connected with a 
rational spline method that allows some curvature between 
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points. This was done to portray the temperature profiles 
of the region mere realistically. 

F. THERMAL FFCFILE TEST AND COMPARISON ANALYSIS 

1 . T J 

The first type of comparison test was conducted with 
rhe criginai rraclc of data used in the generation of the 
predicted thermal profiles. The XET temperature- depth pro- 
files of the Subarctic Water and Transition Zone were com- 
pared with the computed temperature profiles which utilized 
the sea surface temperature of the actual XET trace as the 
only input. Fig. 21 shows the comparison between the 
observed and predicted thermal structures for the northern- 
most XBT station on the track. Appendix F lists the remain- 
ing twenty profile comparisons for the northern regions 
excluding the Subtropic Water. 

The comparisons generally illustrated a strong simi- 
larity in the upper layers. The mixed layer depth along 
with the slight negative temperature gradient in the layer 
were realistically predicted for a majority of cases; inter- 
nal wave motions are assumed responsible for the largest 
variation from the predicted profiles. The upper thermo- 
cline region closely approximated the strong thermal 
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Figure 2 1 . 



Ohserved and Predict ed Thermal P-of-'l-s 
Station 339. 



X3T 
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gradient measured by the XBTs. In the lower therraccline, 
greater variability occurred, assumed due to internal waves 
and thermal advection, resulting in poor comparisons as 
shown in Figs. P.2, F. 1 1 and F. 12 for XBT srarions 34 1 , 350 
and 351. 

The predicted lower regions were not as precise when 
compared to the observed profiles usually due tc the exis- 
tence of a significant temperature inversion. Lack of a 
model for this feature, which was not attempted, clearly 
affected the profiles at depths greater than 130 m. Several 
predicted profiles, however, did attempt tc reproduce this 
structure with the placement of the DT100 variable below the 
BUZ. This was shown in XET stations 343-346 and 348 (Figs. 
F. 4 - F. 7 and F. 9) . 

Between the bottom of the main thermocline and 460 
m, several predicted profiles ccmpletely failed to approxi- 
mate the observed structure. This was caused by the lack of 
variables in the lower layers as well as the computer- gener- 
ated curvature between data points. When the variables 
DT100 and BUZ were very close together, as shown for XET 
stations 348 and 349 (Figs. F.9 and F.10) , the profile would 
bulge with a positive or negative gradient. 
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Otilizing only nine points no define t.he predicted 
thermal profile obviously smooths the fine structure within 
the thermoclins and mixed layer. However, due to the small 
size and transient nature of these features, their removal 
would net severely limit the usefulness of the predicted 
structures . 

In this initial test, since the predicted 
temperature structures were generated and compared with the 
same XBT data from the meridional track, the comparisens 
should be relatively close. The predicted profiles would 
actually represent tha average of all the cbservsd thermal 
structures as the variables reflect the trends along the 
track. 

2. TEST 2 

A more stringent and realistic test was attempted 
that defines the spatial limitations of this thermal pre- 
diction method. Test 2 consists of comparing observed and 
predicted temperature profiles which occurred during the 
same time period as the data generation, but spatially 
removed. A possible indication as to a limit in the effec- 
tive distance would be the observation that the predicted 
profile begins to diverge greatly from the actual 
temperature profile. 
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Th= test XBT data were ottained from Fleet Numerical 



Oceanography Center by a computer search for all available 
XBT drops which occurred in September from 1965 to 1980 in a 
ren degree wide strip of longitude which surrounds rhe orig- 
inal track. Inexplicably, the data received were limited to 
the area south of 44N, excluding the Subarctic Water. Lack 
of tire prevented an attempt to acquire these data whicn 
would be useful, but not entirely necessary, for analysis. 
The temperature-depth data were then selected by noting the 
positions relative to the original data. To provide an 
effective test, a variety of locations in the Transitional 
Zone were chosen. However, these positions were constrained 
by the availability of data in the region and appeared as a 
zonal belt perpendicular to the original track. Fig. 22 
shows rhe locaticn of original and test XBT drops. 

The temperature-depth pcints provided by FNOC were 
then plotted and connected on the computer graphical display 
with the same amcunt of curvature applied to the profile as 
before. Both observed and predicted profiles were plotted 
together again for analysis and comparison purposes. 

Five vertical temperature profiles from the 10-11 
September period were compared with the predicted profiles 
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Figure 22. Lccation cf Original (•) and Tesr 2 (X) X3T 
Drcps. 
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as 



generated from the 



sea surface temperature of 



t ne 



observed profiles. The test produced varied results with 
several features of the computed temperature profiles com- 
paring favorably with the actual thermal structure. The 
largest differences were noted for the two profiles located 
the greatest distance from the meridional track at 145W. 
Figs. 23 and 24 illustrate the comparisons at these loca- 
tions, approximately 400 km east and west from the original 
track of data. The mixed layer depth and upper thermocline 
regions of the predicted profile closely resembled the 
observed profile, but both structures displayed differences 
in the lower regions. 

The eastern most observed profile had a thin, strong 
thermocline which gradually weakened with depth. This con- 
trasted markedly with the strong and sharply-defined thermo- 
cline of the predicted profile. The regions below 150 m 
displayed similar temperature gradients for both profiles. 

For the western most test profile, the only region 
to differ significantly from the predicted structure was the 
lower domain below 130 m. The location of the two close 
variables, BDZ and DT100, along with the computer-gen erated 
curvature forced the predicted profile to diverge from the 
weak temperature inversions. 
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Figure 23. 



Observed and Predicted Thermal Profiles at 
41. 5N, 140. 3W For 101200Z S3? 77. 
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Figure 24 



Observed and 
4 2.3N, 15 0. OW 



Predicted Thermal Pro 
For 110059Z S2F 77. 
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Th€ thrss other test cases, shown in Figs. 25, 26 

and 27, which were located closer to rhe original track, all 
displayed strong similarities between the observed and 
predicted profiles. Except for the fine srructure and pcs- 

( 

sible effects of internal waves cn the depr.h of rhe thermo- 
dins, the computed profiles were very accurate. The 
unusual shape of the observed temperature profile at the 
bottem of the seasonal thermocline, in Fig. 27, and clcse 
approximation by the prscicted profile, illustrates the 
strong ability of this statistical method to produce thermal 

! 

features unique to a certain region. 

To test the spatial limitations of this empirical 
method further, a vertical temperature profile taken at 
Ocean Weather Station "P" during the data-generation time 
period was examined. Since the location of the northern 
most X3T drop cn the track was two degrees of latitude south 
of OWS-P, the results should indicate the reliability of the 
methed outside the test region. Fig. 28 displays the com- 
parison between the observed and predicted prefiles for 
OWS-P. 

Frcm the surface tc the middle of the thermocline, 
the two profiles were very similar, with the mixed layer 
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Figurs 25. 



Ots»rved ’nd Predicted Thermal Profiles a 
41. 9N, 14£.7W For 10 1200Z SS? 77. 
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Figure 26. Otservei and Fredicted Thermal Profiles at 
42. 3N, 144. OW For 102115Z SEP 77. 
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Figure 27 



Observed and Predicted Thermal Pro 
41. IN, 143. 2W For 1 10000Z SEP 77. 



les at 
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Figure 28. Cts«=rved and Predicted Thsrual Profiles at 
fcr 110318Z SEP 77. 



CWS-? 
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depth accurately depicted. However, the lower segmen*: of 
the two profiles differed by 2.0C. With the actual vertical 
temperature prefile being considerably cooler than 
predicted, a different regime must have been encountered. 
Tabata (1961) suggests frem his observations in the region 
that hcritcntal transport cf colder water from the vicinity 
of the Alaskan Gyre, lying northwest of the station, may 
cause a decrease in temperature during the summer months. 

3. TEST 3 

The third test included both spatial and time dif- 
ferences in an attempt tc forecast the vertical thermal 
structure up to three weeks in the future. The only availa- 
ble XBT test da-*-a in the Transitional Region were located 
within three degrees longitude to the west cf the original 
track as shown in Fig. 29. Vertical temperature profiles 
were chosen at weekly intervals and compared with the 
predicted structures. 

At the end cf the first week, the observed thermal 
profile, as shown in Fig. 30, displays a deeper mixed layer 
to a depth of 4C m overlying an extremely strong negative 
temperature gradient. This feature was not represented well 
by the predicted profile as the computed MLD and thermocline 
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Figure 29. 



Location of Original 
Data. 



(•) and Test 3 



(X) XBT 
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variables were shallower in. depth. The shape of the 
observed mixed layer infers strong wind and convective mix- 
ing, so surface weazher charts for the period of 11-18 Sep- 
tember 1977 in the Northeast Pacific ware checked fcr the 
occurrence of a significant atmospheric event. On 17 Sep- 
tember 1977, a cold front associated with a low pressure 
system in the Gulf of Alaska had quickly passed through the 
test region with rain and 35 knot winds. The weather during 
the period of data acquisition, 8-11 September 1977, as 
described before, was mostly cloudy with periods of light 
winds and scattered showers (FNOC, 1977). These contrasting 
atmospheric situations suggest a major limitation in the 
statistical mcdel, as weather effects on the upper vertical 
temperature structure must be prcperly accounted for in the 
predicted profile. 



The lower half of the thermccline on the predicted 
temperature profile closely followed the observed structure 
down to 170 m. Except for the weak temperature inversion at 
200 m, which was not included in the empirical model, the 
lower domain behaved well. 

The next test profile was taker on 24 September 
1977, thirteen days after the data generation for the 
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Figure 30. Observed and Predicted Thernial P’*ofiles ■^o 
182030Z SEP 77 at 43 . 7N , 146. 1 W . 
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predicted structure. Location of this test XBT was on the 
southern hcuncary of the Transitional Zone approximately 200 
ka west of the original track. The predicted temperature 
profile in Fig. 31 clearly resembled the observed vertical 
tharnal structure, but with slightly cooler temperatures. 
Due to the strong winds from the storm a week earlier, the 
observed HID remained deep. However, a transient at 30 m 
depth had formed with a 0.5C increase in sea surface 
temperature as the weather in the region during the second 
week was dominated by a strong high pressure system. The 
observed profile’s thermocline- a Isc deepenad by about 10 m, 
which may be due to internal waves, or it may have been 
caused by wan advection from the south. Immediately below 
the thermocline, the two profiles differed slightly because 
of the fine structure in the region. The predicted 
temperature at 460 m depth remained accurate when compared 
to the observed temperature throughout all the temporal 
tests. 

The final test forecast was conducted for an XBT 
location 209 km to the west of the north- south track on 30 
September 1977, nearly twc weeks after the basic data were 
acquired. Another cold front associated with a low pressure 
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gur€ 31. 



Otssrved and Bredictad Thermal Profiles for 
24C0C0Z SEP 77 at 41. 2N, 147. 5W. 
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system moved through the area on 25-26 September 1977 and 
was immediately followed by a high pressure ridge. The 
storm, with maximum winds of 40 knots, maintained the 40 m 
thick isothermal mixed layer as shown in Fig. 32. The wind 
mixing had increased and deepened the negative temperature 
gradient of the thermocline to a level approximately 10 m 
below the predicted profile. In the lower domain of the 
thermal structure a large difference was noted due to the 
slight temperature inversion of the observed profile and to 
the close location of the two variables, BUZ and DT100, at 
130 II on the predicted profile. The profiles coincided at 
the 460 m depth as determined by the variaole T460. 

The accuracy of the predicted profile after three 
weeks from! the time of data acquisition for the statistical 
analysis shows an overall degradation. However, considering 
that two storm systems passed through the region and the 
summer heating season was approaching its end, the fore- 
casted profile was not excessively poor. 

G. BEVISH CE EHEDICTED TEEBMAL PROFILE RESULTS 

The underlying requirement of the statistical approach 
in predicting the vertical temperature prefile is to acquire 
accurate data that reflects the mean thermal conditions for 
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a particular season and region. Besides the acquisition 

techniques, the importance in precisely defining the 

regional bouncariss is essential. Both temperature and 
saliniry should be used to establish the limits of usable 
data when applying correlation and regression analyses. 
Results of the predicted thermal profiles located outside 
the data generation area, as defined by a certain region, 
displayed a marked degradation in the lower thermocline and 
lower domain regions. The predicted profiles located near 
the area generally portrayed rhe observed temperature pro- 
files realistically within the sparial limits of 455 km 
about the original data. 

The results of the temporal accuracy of the predicted 
profiles indicate a favorable ability to forecast the ther- 
mal structure with certain limitations. The main limiting 
factor was the weather which must be included in the analy- 
sis. Clearly, when high winds pass through a region after 
the vertical temperature data is acquired, the upper thermal 
structure will be altered. Two regions of the thermal pro- 
file that will be affected are the depth of the mixed layer 
and strength of the thermocline. For short periods of about 
one week, this upper domain in the water column will 
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primarily reflect changes due to rhe wind and free convec- 
tive mixing. For time scales cf greater than a week, an 
occasional thermal intrusion from the Subarctic or Subtropic 
Region may be an additicnal factor that could alter the 
thermal structure. Such advective processes would be diffi- 
cult to forecast by the statistical method in regions where 
large meanders and eddies occur, but they should net greatly 
affect a relatively quiescent area such as in the Tran- 
sitional Zens. 

An example cf another problem encountered oy this method 
of thermal structure prediction from the sea surface 
temperature is the process of heating at the ocean's sur- 
face, called the "afterncon effect." Fig. 33 illustrates 
the difference in the predicted profile's structure with the 
observed "afterncon effect." Since the sea surface 

temperature determines the thermal profile in this statisti- 
cal mcdel, an error is imposed in the subsurface structure, 
throughout the entire depth. With the surface heating 

removed, the predicted profile wculd closely resemble the 
observed profile. The "afternoon effect" would greatly 
hinder the ability of remote sensing to determine the sub- 
surface thermal structure by this method. 




Figure 33. 



Example Predicted Profile Error Due to 
"Afternoon Effect". 
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Ihs raethcd of det?r micing mixed layer depth was reviewed 
to ensure that the best possible value was computed and used 
in the empirical scheme. As discussed earlier, the MLD pro- 
vides a base for the variables DT5, DT15, DT50, DTlOO, DT, 
DZ , DTBUZ and DZEUZ to define their values in the predicred 
temperature profile. The ether variables BOT, TBOT, BUZ and 
TBOZ, besides contributing to the construction of the pro- 
file, provida a check to maintain continuity with the MLD- 
determinad variables. The technique used in this thesis to 
compute MLD caused no problems in the determination of the 
thermal profile; however a comparison test of computing 
mixed layer depth by the five different methods, defined 
earlier in Chapter IV, was conducted to ascertain which 
method is the most accurate. 

All methods ware provided with the same SST values from 
the original data of the north-south track and compared with 
tha actual mixed layer depth. The computed depths were than 
correlated to determine the best technique. Results from 
tha test indicated all methods were surprisingly equal. The 
method used was chosen because of the higher correlations 
between variables. Howavsr, the simplist mathod, which was 
the SST-derived mixed layer depth, with a correlation coef- 
ficient of 0.535, could have been used with equal success. 



V. CONCLUSIONS AND EEC C HtN END A TIO NS 
A. CONCLUSIONS 

I Ihs sea surface temperature reflects many of the physi- 

cal processes that govern the oceans. Use of remote sensing 
to map the surface thermal structure of the world's oceans 
accurately and rapidly presently exists, but the determina- 

I 

ti on of precise subsurface vertical temperature structure 
from these surface signatures is unrealized. Possibly only 
in selected ocean regions can insight be gained into deter- 

I mining quantitative relationships for directly inferring the 

( 

subsurface features from the sea surface temperature. In 
I relatively quiet locations with uniformly changing proper- 
ties, such as in the Transitional Zona of the Northeast 

r 

Pacific, many etrcng correlaticns can be observed within the 
, thermal structure. 

The results cf the statistical analysis indicated that 
I the strongest relationships within the vertical temperature 
profile occurred with the thermally-dependent variables 
(TBOT, DI, T£UZ, DTBUZ) as they varied with SST. The 
? depth-related variables (KID, BOT, DZ, BUZ, DZ3UZ) displayed 
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weaker corre laticns with the surface temperature, bur trends 
still could be observed through the thermal structure. The 
use of scatter diagrams to assist in the correlation analy- 
sis cf variables was necessary to establish regional bounda- 
ries and the test fit points for the derivation of 
regression equations. 

An important region within the vertical temperature pro- 
file is the thermocline. In practical terms, its depth, 
strength and gradient strongly influence the propagation cf 
underwater sound. Tne prediction of this structure by use 
of only sea surface temperature would greatly enhance ASW 
capability. Normally, within the Northeast Pacific, the 
correlation between SST and thermocline depth is net very 
high; however, in the Transition Zone with uniformly chang- 
ing structure a reasonable predictability can be obtained. 
The thermocline variables (DT5, DT15, DT50, DT100) all pro- 

vided trends frea which SST-related values were computed to 
define the thermal structure. Overall, the most precisely 
predicted region of the temperature profile was the thermo- 
cline. The depth may have varied due to internal waves, but 
the negative temperature gradient was accurate and realistic 
when compared to the observed structure. 
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In future research cf this mexbod, an expanded ss- cf 
vertical tempeiat’xre data from parallel or bisecting cracks 
should be estatlished as the data base for .statistical anal- 
ysis, The additicn cf an east-west track of data would pro- 
vide latitudinal trends which could also be used in the 
predicted temferature profiles. As in any data aquisiticn, 
a consideration as to the minimum amount of data points to 
determine the thermal structure statistically should be 
made. Meanwhile, a number of no less than 20 profiles over 
a distance of 1000 km with spacing between points not to 
exceed 50 km is recommended. 

ether suggestions for model improvement include adding 
dynamical and thermodynamical variables such as wind speed, 
direction and duration, air temperature, radiation flux, and 
other atmospheric effects which could alter the ocean ther- 
mal structure. Most of these factors can be determined by 
remote sensing; data from other satellite sensors, such as 
an altimeter, SA5 and Coastal Zone Color Scanner which can 
measure sea surface signatures, should also be included. 

To determine the limitation of accuracy for the statis- 
tically-determined profile method, more stringent tests 



should be conducted which increase the spaoial and riaie fac- 
tors. In certain ocean regions one could possibly maintain 
reasonable accuracies within hundreds of kilometers and sev- 
eral weeks duration from nhe acquisition of thermal data. 
Different ccean locations should be analyzed and nested to 
determine the possible regions where nhe thermal profiles 
vary in an orderly fashion which can be easily modeled, 
areas where mescscale features, such as eddies and fronns, 
are consnantly altering the water column by advecnive pro- 
cesses would be poor candidates for this menhod. 

Another re cc crmendation is to test in other seasons of 
the year to determine if the statistical approach would 
apply. Ihe difference in the vertical thermal structure 
from the heating to cooling season is sc great that possibly 
there could be only a certain window of time that allows th^ 
sea surface temperature to define the subsurface structure. 

Finally, a method fcr determining the occurrence of 
large temperature inversions in the Subarctic and Transition 
Regions is needed. This important thermal feature, because 
of its varlaticn in magnitude, depth and location, affects 
greatly the propagation of sound by forming a sound channel 
immediately belcw the the r race lin s. If this structure can be 
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predicted then the full potential of the statistical niethcd 
of determining subsurface thermal structure from sea surface 
temperature would be realized. 
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SFPENDIX A 



THERMAL VARIABLES VERSUS LATITUDE/DISTANCE ALONG 
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Figure A. 1. 



SST Versus 



Latitude/Distance. 
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MIXED LAYER DEPTH (M) 
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Figure A. 2. 



MLE Versus Latitud e/Distarce. 
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Figure A. 3. DI5 Versus Latitude/Disr anca. 
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Figure A. 4. 



CT15 Versus Laritude/Disrance. 
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Figure A. 5 
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Figure A. 6. 



EH00 Versus La tit ude/Diszar ca . 
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Figure A. 8. 
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Figure A . 1 1 . 



BUZ Versus Latitude/Disrancs. 
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Figur-s A. 12. TEOZ Versus Latit ude/D isrance. 
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AEPENDIX B 



BICKEDICAl COMFUTEE PROGRAilS (DIXON AND BROWN 

1979) 



Ths EIC.'iED statistical computer programs wars utilizsi 
because of the simple and flexible programming it provided 
in analyzing the thermal data. In constructing the vertical 
temperature profiles, the statistical analysis required cor- 
relation and regression equation computation and scatter 
diagram plotting of the thermal variables. 

The program initially selected was 'S-.epwise Regression 
P2R.' This program ccirputes simple statistics, i.e., 
mean, standard deviation, etc., correlation coefficients and 
linear regression equations for any number of variables. 
The statistics used in the computations are defined as: 



X = ^Xj/N 





* the correlation coefficient between two variables 



r - 
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In simple linear regression 



analysis , 



where there 



only one independent variable, the coefficients are deter- 
mined by a least squares method. The form of the equations 
ar e: 

Y = A + Bx (where A and B are the linear coefficients) 

The EICMID computer program used to print the regional 
scatter diagrams is titled ’Bivariate (Scatter) Plots 
P6D.’ Recressior. equations for regressing variable x on y 
and y on X can be requested along with the intersections of 
the regression lines with the axis of the plot. This was 
useful in analyzing the data to determine the best fit for 
the various regions. 
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Figure C.1. Begional Scatter Diagram - SST Versus MLD. 
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Figure C.8. 
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SST Versus DT. 
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Figure F . 2. 



Observed and Predicted Thermal Profiles at XET 
Sta. 341. 



169 



Figure 




F. 2. 



Observed and Predicted 
Sta. 342. 



Thermal Profiles a 



t XET 



170 



DEPTH (METERS) 

500 -4i0 ‘400 -350 -300 -350 -200 




Figure F.4. 



Observed and Predicted 
Sta. 343. 



Thermal 



Profiles a 



t XBT 



171 



DEPTH (METERS) 

SOO -490 -400 -390 -300 -290 -200 




Figure F.5. 



Observed and Predicted Thermal 



Profiles at 



XBT 



172 



DEPTH (METERS) 

aoo ~<30 -^00 -350 -300 -250 -200 




Figure F.6. 



Observed and Predicted Thermal 
£ta. 345. 



Profiles 



at 



X3T 



173 




I 





Figure F.7. 



Observed and Predicted Thermal Profiles at XBT 
Sta, 346. 



174 







/ 




o 



Fi gure 



o 

«- 



O ' 
O- 



O ' 
*r>- 




O- 

c/> I 
cx 

ui _ 

5 w 



CL 

O o_ 
« 

I 



o 

W- 



o 

o. 



Sj 



LEGSND 
a 02Si.WED 



o PREDiaED 



O ' 

I ^ 



2.5 5 7.5 10 12.5 15 17,5 

TEMPERATURE CDEGSEES O 



20 



22J 25 



F . 6. 



Observed and Predicted Thermal Profiles an XET 
Sta. 3U7. 



175 




\ 




DEPTH (METERS) 

SOO -M90 ’400 -3S0 -900 -290 -200 




Figure F. 9. 



Observed and Predicted Thermal Profiles at XBT 



176 



DEPTH (METERS) 

300 -430 -400 -330 -300 -230 -200 




Figure F . 10 . 



observed 
XET Sta. 



and Predicted Thermal Prcf 
349. 



s at 



177 




Figure F. 11 . 



Obssr ved 
XEI Sta. 



ard Predicted Thermal 
350 . 



Prof il es 



at 



178 



DEPTH (METERS) 

-900 -430 -400 -330 -300 -330 -300 




Figure F . 12 . 



Observed 
XET Sta. 



and Predicred Thermal Profiles 
351. 



at 



179 




gurs F . 13 . 



Observed and Predicted Thsraal P’-of-'le 
XET Sta. 352. 



180 




Figure 



F . 14. 



0 bs9"^ved 
XEI Sta. 



and Predicted Thermal Profiles 
353. 



at 



181 




Fi gu r € F . 15 . 



Obc9rv5d 
XET Sta. 



and Predicted Thermal Profiles 
354. 



182 




f 






i 





Figure F . 16 . 



Cbssrved and Predicted Thermal Profiles at 
XBT Sta. 355. 



183 




Figure F . 17 . 



C bssr V ed 
XET sta. 



and Predicted Thermal Profiles 
356. 



af 



184 




gurs 



F. 18. 



otser v-3d 
XET Sta. 



and Predicted Thsrmal 
357. 



Prof lies 



an 



185 




Figure F . 19 , 



Observed 
XET St a. 



and Predicted Thermal Profiles at 
358 . 



186 




Figure F . 20 . 



Ctssrved 
XET Sta. 



and Predicted Thsrinal Profiles ar 
359 . 



1 87 



I 




LIST OF REFERENCES 



Bath-sn, K.H. , "Or. ths Ssascnal Changes in nhe Depth cf the 
Mixed Layer an the North Pacific Ocean", Journal of 
Ge ophysical Research, v.77, n.36 , 7 138-7 15J7~tJ'?cenDer , 1972. 



Ba th ytker mcgraph Log Sheets 113-122 for SILA.S BENT cruise 
343722, Phase III, September 1S77. 



Bernstein, H., Breaker, L.^ and Whritner, R. , "Califcrnia 
Current Eddy Formation; Shin, Air and Satellite Results", 
Science, v.195, 353-359, January 1977. 



Clancy, R.K., and Martin, P.J., "Synoptic Forecastina of the - 
Oceanic Mixed Laver Using ths Navy’s Operational 
Environmental Data Base: Present Capabilities and Future 

Applications", Bulletin of the American Metecrc Icqica 1 
Societi, v.62^ n7^~77TJ-7^47-Jun~T9HlT 



Cox, S.A., Satellite Ap plication to Acoust ic Pr edi ction 
Sv stems, Master‘'’sT1i9sis7~T3av^ Postgraduate School, ~ 
dontirey, 1982. 



Denman, K.L., and Miyake, M., "Upper Layer Modification at 
Ocean Station Para; Observations and Simulations", Journal 
of Physical Ocear. cqr a phy , v.3, n.2, 185-196, April T773T 



Dixon, W.J., and Brown, M.E., Biomedical Computer Proarams 
P- Series, EMDP-79, University 5r raXiTornia ?ress, ~5erTceTey , 
CaIiIcfnia7~l-ddU, 1979. 



Dodimead, A.J., Favorite, F., and Hi 
Oceanography of the Subarctic Pacific 
lit, ■Paper“o recafld for‘"t7e”Tnt ef nati 
Fisneries Commission, 1 October 1962. 



rano, T. , Rev iew 
Reqi on , V . 7, f i 
onaX dorth Pacif 



gure 

1C 



Fleet Numerical Oceanography Center, Sea Surface Temperature 
Evaluation for 1C1200Z September 1977. 



Fleet Numerical Oceanography Canter, Surface Weather 
Analysis for September 1977 



Frankigncul C., "Low Frequency Temperature Fluctuations Off 
Bermuda", Journal of Geophysi cal Research, v.86, n.7, 

65 22-652 8,-7uIy~79STT. ^ 



Kenycn, K.E., "A Shallow Northeastern Current in the North 
Pacific", Journal of Gsoohysical Research, v.86, n.7, 
6529-6536 , “JuIy“79HT. 



Kenycn, K.E., "The Surface layer of the Eastern North 
Pacific in Winter", Journal or Geophysical R ese arc h, v.83, 
n.12, 6115-6122, Dece¥5ef "797B. 



138 



I 




Lsgeckis, R., "A Survsy cf Worldwide Sea Surface Temoerature 
Fronts Detscosd by Environmental Satellites", Jour na 1 cf 
Geo£hysical Rese arc h , v.83(C9), 4501-4522, Seprsmoer 



Lageckis, E. , and Gorden, A.L., "Satellite Observations cf 
The Erazil and Falkland Currents - 1975 to 1976 and 1978", 
De ep -S ea Research, v. 29, 392-405, March 1 982. 



Miller, I., and Freund, J.E., Pr obab il iry and Star isr ic s f or 
Engineers, Prenr ice-Hall, Inc., ‘Englawood Clifis, TIiw ~ 
3ersey7~35 1-328, 1977 . 



Monin, A.S., and Fedorov, K.N., "The Fine Srrucrure of the 
Upper Lavers cf rhe Ocean", Atmospheric and Oceanic Physics, 
V.9, n.4, 442-444, July 19727 



Renner, J.A., "Fishina I nf crmaticn" , Southwesr Fisheries 
Center - La Jclla, California, n.9, 1-13, September 1977. 



Roden, G.I., "Aspects of the Mid-Pacific Transition Zone", 
Journal of Gecphysical Research, v.75, n. 6 , 1097-1109, 

Pe of uary“3 9707 



Roden, G.I., "Cn North Pacific Temperature, Salinity, Sound 
Velocity and Density Fronts and Their Relation to Wind and 
Energy Flux Fields", Journal of Physical Oceano era ohy , v.5, 
n. 4,''557-571 , CctoberT9757 



Roden, G.I., "Shallow Temperature Inversions in the Pacific 
Ocean", Journal cf Gecphysical Research, v.69, n.14, 

28 99-29 117"3uly T56ir7 



Tabata, S., "Characteristics of Water and Variations cf 
Salinity, Temperature, and Dissolved Oxygen Content of the 
Water at Ocean Weather Station ' P' in tne Northeast Pacific 
Ocean", Journal cf the Fisheries Research Board of Canada, 
V. 17 (3) ,“53=j707“17^-a. 



Tabata, S., "Comparison of Observations of Sea Surface 
Temperatures at Ocean Station ? and NOAA Bouy Stations and 
Those Made by Merchant Ships Traveling in Tneir Vicinities, 
in the Northeast Pacific Ccean", Journal cf A pplie d 
^tecrolo_gy, v. 17, 374-385, 1978 . 



Tabata, S., "Variability 
Ocean Station * 

Transactions _ 

3H7-TnH7~3 un e~l 9F57 



cf Oceanographic Conditions at 
•P' in the Northeast Pacific Ocean", 

<jt the Royal Socie t y cf C^ada, v.3. Series 4, 



Tabata, s., Ecsten, N.E.J., and Eoyes, F.M., "The Relation 
Between Wind Speed and summer Isothermal surface Layer of 
Water at Ccean Station P in the Eastern Subarctic Pacific 
Ocean". Journal cf Gecphysical Rese arc h. v.70, n.16, 
3867-38777~SugXsf”17^57 



189 




i 

i 



i 




Tabata, S., and Kimber, P.M., "Satellina 0 bssr var:icns of Saa 
Surface T smpsr ature Patterns Off the Pacific Coast of 
Canada'*, Pacific Marine S c ien ce Reoort 79-19, Institute of 
Ocean Sciences7 Patrrcia~'Eay,“BrltasH~C^uir5aa, Canada, 

1-77, 1979. 



Tally, J.P., "Oceanographic Regions and Assessment of 
Temperature Structure in the Seasonal Zone of the North 
Pacific Ocean", Jcurnal of the Fisheries Research Board of 
^nada, v . 2 1 (5) ,"'95'T-‘97 07“Ssptem'Eer^'9&‘r. ~ 



Tally, J.P., and Giovando, L. F. , "Seasonal Temperature 
Structure in th® Eastern Subarctic Pacific Ocean". Marin e 
Distribut icns , The Royal Society of Canada Special 
FuIIica^Icn'llc. 5, 1 0-35, 1963. 



Uda, Michitaks, "Cceanoarathy of the Subarctic Pacific 
Ocean", Jcurnal of the Fisheries Research Board of Canada, 

V. 20 (1) ,"11^=1797~19^3. 



White, W. , Bernstein, R. , McNally, G. , and Pazan, S. , "The 
Thsrmoclins Response to Transient Atmospheric Forcing in the 
Interior Mid-Latitude North Pacific 1976-1973", Journal of 
Ph ys i ca l Cceanoar^hv, v. 10, 372-384, March 1930. 



190 



INITIAL DISTRIBUTICN LIST 



No. 



1. Dafanss Technical Ir.fcrmaticn Cenxar 
Cameicn Station 
Alexandria, Virginia 22314 



2. Library, Cods 0142 

Naval Postgraduate School 
Monterey, California 93940 



3. Chairman (Code 68Mr) 

Department of Oceanography 
Naval Postqraduate School 
Monterey, CA 93940 



4. Chairman (Code 6 3Rd) 

Department cf Meteorclcgy 
Naval Postgraduate Schcci 
Monterey, CA 93940 



5. Superintendent 
Attn: Code 68Jg 
Naval Postqraduate School 
Monterey, California S3940 



6. Superin terdent 

Attn: Code 68Du 
Naval Postgraduate School 
Monterey, California 93940 



7. Director 

Naval Oceanography Division 
Naval Observatorv 
24th and Massachusetts Avenue NW 
Washington, C. C. 20390 



8. Commander 

Naval Slectrcnics Systems Command 
Attn: PME-124-60 

Naval Electronics Systems Command Headquarters 
Washington E.C. 203oC 



9. Commander 

Cceancaraphic Systems, Pacific 
Eox 1390 

Pearl Harbor, HI 96360 



10. Commander 

Attn: Cede 531 

Naval Ocean Systems Center 

San Diego, C A 92152 



1 91 



copi 

2 

2 

1 

1 

2 

2 

1 

1 

1 

1 



11 



1 



. Chief cf Naval Oce rations 
Attn: OE951F 
Navy Cecartient 
Wasningtcn C.C. 20350 



12. Commanding Officer 1 

Fleet Numerical Oceanography Center 

Monterey, CA 93940 

13. Commanding Officer 1 

Naval Ocean Research and Develooment Activity 
NSTL Station 

Bay St. Lcuis, MS 39522 

14. Commanding Officer 1 

Naval Environmental Prediction Research Facility 
Monterey, CA 93940 

15. Chainran, Oceanography Department 1 

U.S. Naval Academy 

Annapolis, MD 21402 

16. Chief of Naval Research 1 

800 N. Quincy Street 

Arlington, VA 22217 

17. Office of Naval Research (Cede 480) 1 

Naval Ocean Research and Development Activity 
NSTL Station 

Bay St. Lcuis, MS 39522 

18. Scientific Liason Office 1 

Office of Naval Research 

Scripps Institute of Cceano graohy 
La Jclla, CA 92037 ‘ 

19. Library 1 

Scripps Institute of Oceanography 

P.O. Box 2367 
La Jclla, CA 92037 

20. Library 1 

Department cf Oceanography 

University cf Washington 
Seattle, KA 98105 

21 . Library 1 

CICESE 

P.O. Box 4803 

San Ysidre, CA 92073 

22. library 1 

School cf Oceanograohy 

Cregen State University 
Corvallis, OR 97331 



192 



? 



23. LT T.A. Hcw€ll 1 

Naval Wastern Ocsano graohy Center 
Box 113 ^ i 

Eearl Harhor, HI 96860 



24. Comniander 1 

Naval Oc5ancgraphy Ccmmand 
NSTL Station 
Bay St. Louis, MS 39522 



25. Commanding Officer 1 

Attn: Dr. William Jotst, Code 7300 
Naval Oceanographic Office 
NSTL Station 
Bay St. Loius, MS 39522 



26. Commanding Officer 1 

Attn: Cods 510 1 
Naval Research Laboratory 
Washingron D.C. 20375 



193 




! 

I 



1 





201550 



Thesis 
H8247 
c. 1 



""^A^statistical approach 
for determining subsur- rc 
face thermal structure -- 
from sea surface temp- - 
erature in the nort 




H8247 i^°'"®5^^statistical approach - 
*'■ for determining subsur- 

face thermal structure 
from sea surface temp 
erature in the north- 
east Pacific Ocean. 



